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We report on the status of our kaon semileptonic form factor calculations using the highly- 
improved staggered quark (HISQ) formulation to simulate the valence fermions. We present 
results for the form factor /+"(()) on the asqtad Nf = 2+1 MILC configurations, discuss the 
chiral-continuum extrapolation, and give a preliminary estimate of the total error. We also present 
a more preliminary set of results for the same form factor but with the sea quarks also simulated 
with the HISQ action; these results include data at the physical light quark masses. The improve- 
ments that we expect to achieve with the use of the HISQ configurations and simulations at the 
physical quark masses are briefly discussed. 
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1. Introduction and methodology 

The study of exclusive semileptonic decays of D and K mesons provides a way of extracting 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements |V C( ^ CJ )| and |V„. S | with errors competi- 
tive with those obtained using other methods such as leptonic decays, neutrino-antineutrino inter- 
actions, or % decays. Comparison of the values obtained with different methods could reveal new 
physics (NP) effects, and comparison of the shape of the form factors describing those exclusive 
decays with experiment can provide a check of the lattice methodology employed. 

Our program includes analyzing K — > nlv and D — > K{n)lv semileptonic decays at zero as 
well as non-zero momentum transfer. In these proceedings we focus on the status of the K — > %lv 
analysis at zero momentum transfer, which one can combine with experimental data to extract the 
CKM matrix element \Vus\- The limiting error currently comes from the lattice determination of 



the form factors defined in (1.1) [ffp. A precise determination of \V US \ provides stringent tests of 
first-row unitarity and gives information about the scale of new physics [Q]. 

The theory input needed to extract the CKM matrix elements from exclusive semileptonic 
widths are form factors parametrizing the corresponding hadronic matrix elements: 

{P2\V^)=f^(q 2 )(p Pl +pr 2 -Ay+f^(q 2 )A^ (1.1) 

where = (m 2 Pi — m 2 p 2 )q^ / q 2 , q = p Pl — pp 2 , and V is the appropriate flavor-changing vector 
current. We obtain the needed form factor /f 71 (0) using the relation /jf 71 (q 2 ) = m f m \ (n\S\K)(q 2 ), 

m K ™n 

and the fact that /+(0) = /o(0) due to the kinematic constraint. This method [J3|] allows us to 
eliminate the need for a renormalization factor and to extract the form factor from three-point 
correlation functions with insertion of a scalar current instead of a vector current. 

The momentum transfer of the three-point functions is tuned to zero or very close to zero 
using twisted boundary conditions to inject external momentum The general structure of the 
correlation functions is given in Fig. [l[ We consider either a moving % (d = d\ = and 6 2 ^ 0) or 
a moving K (6 = d 2 = and B\ / 0), for K ->• nlv. 



2. HISQ valence fermions and asqtad Nf = 2 + 1 configurations 

For the first analysis we use the asqtad Nf = 2 + 1 MILC configurations and the HISQ action 
to simulate the valence quarks. The strange valence-quark mass is fixed to its physical value and 
the light valence-quark mass is fixed so "l'h y f HISQ ^ = '"'i lv f asqtad ^ . The parameters of the ensembles 

m s ' (HISQ) m s (asqtad) 

included in this calculation and the details of the simulations we perform are collected in Table [l|. 

In order to extract the value of the form factor f+ K (Q) = /^(O), we perform a simultaneous fit 
of the relevant three- and two-point functions. The correlator fits and the multiple checks performed 
on their stability under the change of parameters, time range, number of exponentials in the fitting 
functions, and correlators included, was described in last year's conference proceedings [§]. The 

* Speaker. 

"'"Operated by Brookhaven Science Associates, LLC, under Contract No. DE-AC02-98CH10886 with the United 
States Department of Energy. 

^Operated by Fermi Research Alliance, LLC, under Contract No. DE-AC02-07CH1 1359 with the United States 
Department of Energy. 
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Figure 1: Structure of the 3-point functions 
needed to calculate fft K ^ [/jf K \. Light quark 
propagators are generated at fsource with ran- 
dom wall sources. An extended charm [strange] 
propagator is generated at T + fsource- 



~ a (fm) 


ami /anih 


Nconf 


N sources 


N T 


0.12 


0.020/0.050 


2052 


4 


5 




0.010/0.050 


2243 


4 


8 




0.07/0.050 


2109 


4 


5 




0.05/0.050 


2098 


8 


5 


0.09 


0.0124/0.031 


1996 


4 


5 




0.0062/0.031 


1946 


4 


5 



Table 1: Aqstad nsembles and simulation details, arrih is 
the nominal strange-quark mass in the sea sector, N SOU rces 
is the number of time sources, and Nj the number of sink- 
source separations. 



conclusion is that we find it very difficult to make changes in the fitting procedure that change the 
fit results outside the one sigma range. The only change with respect to last year's fits is that in our 
final correlator fits we are including functions with the momentum injected in both the n and in the 
K. The results from these combined fits for the different ensembles are shown in Fig. ||[ Statistical 
errors are very small, ~ 0.1 — 0.15%. 

2.1 Chiral-continuum extrapolations 

We need to extrapolate our results to the continuum limit and the physical light-quark masses, 
and also adjust for the mistuning of the strange sea-quark mass in the asqtad Nf ■ = 2 + 1 MILC con- 
figurations. Our plan is to perform these extrapolations using partially quenched staggered chiral 
perturbation theory (S#PT) at NLO plus regular continuum ^PT at NNLO. Addressing staggered 
effects at NLO should be enough to achieve the sub-percent precision that we target. 

In the continuum, the form factor is given by ff x (0) = 1 + fi + + + . . . , where, according 
to the Ademollo-Gatto (AG) theorem, the chiral corrections f%, fy, f$ ... go to zero in the SU (3) 
limit as (nr K — m 2 n ) 2 . This means that at NLO there are no free low-energy constants and fi is fixed 
in terms of experimental quantities. At finite lattice spacing, however, we would have violations of 
the AG theorem due to discretization effects in the continuum dispersion relation needed to derive 
the relation between f^ n (0) and the correlation functions we are generating. 

The general structure of the fitting function we plan to use for our chiral fits is thus 



/f (0) = l+f 2 Qsta8 \a) + Cf> ( f ) (logs) 



a 



7i 

+4 (ml -m 2 K ) 2 



(2.1) 



(a) A") 



'(1) 



and Cg 1J °= C12 + C34 — are free parameters to be fixed by the 



where the constants C\ , C L 
chiral fits. The L' t s are the usual 0(p 4 ) low-energy constants (LEC's), and the Cy are ^(/? 6 ) LEC's 
denned in The function /^'^'(a) is the NLO partially quenched S#PT expression, which 
incorporates the dominant lattice artifacts from taste breaking. To a very good approximation, there 
are no free parameters in that function. The taste-splitting and the taste- violating hairpin parameters 
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are already available for both the valence and the sea quarks from lattice calculations for asqtad 
or HISQ fermions, respectively. We are also in the process of calculating the taste-splittings for 
the mixed mesons, those made of one sea fermion and one valence fermion. The only parameter 
that we do not get from other calculations is an extra taste-violating haiipin parameter that appears 
due to the fact that we have a mixed action. We need to leave this quantity as a free parameter of 
the fit, although we expect its impact to be small. We take the NNLO contribution, fa, from the 
calculation in Finally, we try to include terms proportional to a 2 with free parameters Q (to 
take into account the violations of the AG theorem at finite lattice spacing) and Cg (to account for 
the residual a 2 dependence at higher orders). 

Since we have not completely checked the NLO partially quenched S^PT calculation, here 
we use a simplified version without hairpin terms. We also approximate the taste-splittings of the 
mixed action mesons by the average of the sea and the valence values, A m j x = (A sea (asqtad) + 
A va i ence (HISQ))/2 (we checked that using the preliminary values for the correct splitting does not 
change the extrapolated ff n by more than 0.1%). With these simplifications, we tried several 
variations of the fitting function in (2A): fixing the LEC's L,'s to their value from the global fit 
in [§J, fixing them to the value from the fits in or leaving the L,'s as free parameters in the fit 
with prior central values equal to the results in [Jsj] and varying the prior widths from the errors in [||] 
to an order of magnitude larger; including only the term proportional to Q , the one proportional to 
Q , or both, etc. The extrapolated value for ff K (0) has statistical errors between 0.2% and 0.3% 
in all cases and the different results agree with each other within one statistical a. The violations 
of the AG theorem are around 0.32 - 0.15% for a 0.12 fm and 0.15 - 0.1% for a « 0.09 fm. A 
typical example of the fits we have performed is shown in the left side of Fig. |2| 

In order to check the impact of the choice of fitting function in the extrapolation, we have 
also done a number of fits replacing the NNLO continuum %?T functions by a NNLO analytical 
parametrization 

/p(o) = i+ f 2 Q ^- +c<f (^y+if^-^^CW+^CW 

+ C%\nm n f\n{ml/ l l 2 ) + C%\nm n f + C l g ) 



(2.2) 



where f^ stag - is the same partially quenched NLO S#PT expression as in (2.1) and the Cf , df\ 
and Cg with i= 1 — 4 are the free parameters in the fit. Again, we tried variations of the functional 
form in (^2f), turning on and off different NNLO and a 2 terms, and parametrizing them in different 
ways. All the fitting functions we tried in this category gave results within one statistical a of each 
other. In the right-hand side of Fig. ^| we plot an example of these fits. Both methods to describe 
NNLO contributions give results that, again, agree within one statistical a. 



3. Preliminary results with HISQ valence quarks and HISQ Nf = 2 + 1 + 1 
configurations 

The second stage of our semileptonic decay program is the study of decays with q 2 = for 
D and K mesons using HISQ for the valence and the sea quarks, i.e., simulating on the HISQ 
Nf = 2 + 1 + 1 MILC configurations. 
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Figure 2: Left-hand plot: example of chiral-continuum extrapolations using the fit function in Eq. (2.1) 
with Cg = and free L,'s and priors equal to the results in Ref. [||]. Right-hand side: example using the fit 
function in Eq. (2.2) with Cg = = = 0. Errors are statistical only, from 500 bootstrap ensembles. 



The setup and methodology of the calculation is common to the asqtad Nf = 2 + 1 calculation 
and is described in Sec. [j]. Aside from the reduction in discretization errors provided by having 
HISQ sea quarks, the main improvement of this calculation respect to the asqtad Nf = 2+1 one is 
that we will include ensembles with physical light-quark masses. The parameters of the ensembles 
we plan to include in this analysis, as well as the status of the runs, are shown in Table ||. The plan 
is to have around 1000 configurations per ensemble and 4-5 different source-sink separations per 
ensemble at three different values of the lattice spacing. For the key ensembles we are generating 
data for 8 time sources, and, for the remaining ensembles, 4 time sources. For the light and strange 
valence-quark masses we will use the physical values. For the charm-quark masses, we will simu- 
late at a value equal to the sea charm-quark mass in addition to the current estimate of the physical 
one, to allow for later corrections of am% hys . 



w a (fm) 


ami /anif, 


am c 


Volume 


N CO nf available 


N sources 


N T 


% run completed 


0.15 


0.035 


0.831 


32 3 x 48 


1020 


8 


5 


100 


0.12 


0.200 


0.635 


24 3 x 64 


1053 


4 


4 


100 




0.100 


0.628 


32 3 x 64 


1020 


4 


4 







0.035 


0.628 


48 3 x 64 


460 


8 


4 


50 


0.09 


0.200 


0.440 


32 3 x 96 


1011 


4 


4 







0.100 


0.430 


48 3 x 96 


1000 


4 


4 







0.035 


0.432 


64 3 x 96 


497 


8 


4 






Table 2: HISQ ensembles and simulation details. N SOU rces is the number of time sources, and Nj the number 
of sink-source separations for which we have generated data. The number of configurations, N con f available 
and the status of the runs (last column) correspond to July 2012. 



In these runs we generate the correlation functions needed for the calculation of the form 
factor at zero momentum transfer for both K — > nlv and D — > K(n)lv. The energies of the pions 



5 



K semileptonic decay form factors with HISQ valence quarks 



E. Gamiz 



^ ^ 



0.12 fm; 111=0.2111^ 
0.15 fm physical n 
0.12 fm physical n 
0.12 fm; m=0.1in 




continuum NLO 

continuum NNLO (fit. asqtad data only) 
a = 0. 1 2 fm (N f = 2+ 1 asqtad configurations) 
a = 0.09 fm (N - 2+1 asqtad configurations) 
a = 0.12 fm (N f = 2+1+1 HISQ configurations) 
a = 0. 1 5 fm (N f =2+1+1 HISQ configurations) 
Extrapolated value from NNLO asqtad fit 



ir, Pp r 



(r m ) 



Figure 3: Deviation of our data from the continuum 

dispersion relation prediction. Figure 4: Form factor f+ K {0) as a function of the 

K mass from the asqtad Nf = 2 + 1 and the HISQ 
Nf = 2+ 1 + 1 calculations, together with the results 
from a fit to asqtad Nf = 2 + 1 alone, also plotted on 
the left-hand side of Fig. 2. 

and kaons generated on those ensembles show very little deviation from the continuum dispersion 
relation, see Fig. ||. The points with larger errors in that plot correspond to the energies needed 
to inject momentum in a K to get q = in K -+ nlv when we have physical light-quark masses. 
Thus, for physical light-quark masses, moving pions will give us significantly smaller statistical 
errors than moving kaons in K — > nlv decays. 

We fit the correlation functions generated on the HISQ configurations for the K -+ nlv de- 
cays using the same fitting functions and following the same strategy described in Sec. |]. The 
preliminary results from these fits are shown in Fig. f|, together with the data generated on the 
asqtad configurations and the results from the fit to the asqtad data alone (left plot in Fig. ^). The 
Nf = 2+1 + \ HISQ data are very close to the continuum line obtained by fitting the Nf ■ = 2 + 1 
asqtad data. In particular, the point corresponding to the ensemble with physical quark masses and 
a f» 0.12 fm lies right on top of the extrapolated value we got from the asqtad fit. This seems to in- 
dicate that the discretization effects in the HISQ data are going to be smaller than in the asqtad data, 
as expected. The statistical error of the physical mass point is ~ 0.2%, larger than the 0. 1 — 0. 15% 
error we got for larger masses, but of the same order as the extrapolated value. 



4. Conclusions 

We have nearly completed the calculation of f+ K (0) at two different values of the lattice spac- 
ing using the asqtad Nf = 2+ 1 MILC configurations. The last step towards finishing the calcula- 
tion is checking the NLO partially quenched S^PT expressions and completing the error budget. 
We estimate that the total error is going to be between 0.35 — 0.5%, dominated by the statistical 
and extrapolation errors (0.2 — 0.3%) and the uncertainty associated with the deviation of amf 
from the physical value (~ 0.2%). We are also investigating the impact of subleading errors, such 
as finite volume effects. The total error will be competitive with current state-of-the-art calcula- 



tions [10, 11]. 
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The dominant two errors in the calculation on the asqtad Nf ■ = 2 + 1 configurations will be 
reduced in the next step of our program namely, the calculation on the HISQ Nj ■ = 2 + 1 + 1 
configurations, for which we have shown preliminary results here. Having data at the physical 
quark masses will reduce the statistical and extrapolation errors, as well as the one associated with 
the choice of chiral fitting function. Discretization errors are also considerably smaller for the HISQ 
action than for the asqtad action, as explicitly seen in Fig. ||. Finally, the strange sea-quark masses 
are much better tuned on the HISQ ensembles, and we are including the effects of the dynamical 
charm-quark. 
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